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Abstract 

There are two formulations at non-zero chemical potential; one is the formulation that a Lagrangian 
includes a chemical potential, the other is the formulation that a Lagrangian does not include a chemical 
potential. The existence of two formulations makes a calculation complicated. The results from those 
formulations are not corresponding directly. This discrepancy exists in the imaginary time formalism and the 
real time formalism. However, since this is essentially caused by a difference of a frequency, the discrepancy 
vanishes by modifying the Fourier transform. We show a calculational procedure with a spectral function 
to understand this. 



1 Introduction 

As shown in Ref. [1], there is a cumbersome puzzle at finite temperature and chemical potential. The puzzle is 
caused by the existence of two formulations in the real time formalism at non-zero chemical potential[T][2]. This 
raises a problem on a specific calculation. For example, a problem appears on a calculation of a summation at 
finite temperature and chemical potential. However, if the imaginary time formalism [3] [1] is consistently used, 
this problem does not emerge. It emerges when a method with real time is adopted. In particular, it emerges 
more directly in the real time formalism [Ij [2]. For this reason, clarifying details of the puzzle is important to 
maintain the consistency between imaginary time and real time. 

There are two formulations at non-zero chemical potential; one is the formulation (A) that a Lagrangian 
includes a chemical potential i5 , the other is the formulation (B) that a Lagrangian does not include a chemical 
potential[5]. A problem does not emerge as long as only either is used in the real time formalism. However, 
when one verifies whether a result agrees with a one from the imaginary time formalism, a problem emerges. 

Since the formulation (A) is used in the imaginary time formalism, the formulation (A) with real time agrees 
with a result from the imaginary time formalism straightforwardly. On the other hand, the formulation (B) 
with real time does not agree with a result from the imaginary time formalism directly. Thus, most results in 
the formulation (A) and (B) are not corresponding. 

We show topics that relate to a spectral function Since it is useful to use a spectral function at finite 
temperature and density's" fO" fTO^, an attention for the use might be important. We show formulations (A) and 
(B) in section [21 In section [31 we explain the cause of the discrepancy intuitively. In section [31 we perform a 
specific calculation using the formulation (A) and (B), and we show the difference on a calculational procedure. 



2 Fermion propagator by a spectral function 
2.1 Formulation (A) 

The ensemble average of an operator at temperature T = 1//3 and chemical potential /x is defined as 

{6)p = Z-itr[e-^(^-^*)0], Z = tr[e-^(^-^*)], (1) 
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where H and N are a Hamiltonian and a number operator. Since we use a fermion field and the time contour 
—oo ^ +00 (in our purpose it is sufficient to consider only this time contour), the time ordered propagator is 
defined as 



S{x,y) ^9{xo - yo){ip{x)tp{y))i3 - 9{yo - xa){ip{y)ip{x)) ,3 

=9ixo - yo)S> [x, y) + e{yo - xo)S< [x, y), (2) 

The advanced and the retarded green functions are defined as 

SA{x,y) = ie{xo - yQ){S>{x,y) - S<{x,y)), (3a) 
SR{x,y) = -i9{yo - xo){S>{x,y) - S<{x,y)). (3b) 
In addition, we define the thermal green function as 

SfsiT, X, r', y) = 9iT - t')(^(t, a;)^(T', y))^ - 9{t' - t)(^(t', y)^(r, x))p (3c) 

where r is imaginary time (t = it, < r < ^). 

We consider a case that the time evolution is given by 

-^{t, x) = e'(^-''*)*V;(t = 0, a;)e-*(^-^*)*. 
By this definition, the Kubo-Martin-Schwinger (KMS) condition for a c-number fermion field is 

m = -^{t - ip). (4) 

The KMS condition for 5'>(i, a;) is 

S>{t,x)^~S<{t + il3,x), (5) 

and S^{p) has the relation, 

S>{p)^ I dt f d^xe'P''S>{t,x) = -e^P°S<{p). (6) 

'J —00 'J 

We define the spectral function in the momentum representation as [5] [7] 

p{p)^S>{p)^S<{p). (7) 

Using this spectral function and dH]), 

S>ip) = (1 - nF{po))pip), S<ipo) = -nFipo)p{p), (8) 

with 



The advanced, the retarded and the thermal green functions can be expressed by a spectral function. An explicit 
expression for spectral function is given by [3 



p(po) - 27rZ-i ^ e-'5^"(l + e^(^"-^"))(m|V(0)|n)(n|7^(0)|m),S(po + - K^), (9) 



where Kn — En — fJ,Nn- En and Nn are eigenvalues of the operator pQ and iV„ for eigenstates |n). We ignore 
spatial components for simplicity. Using this expression, the advanced, the retarded, and the thermal green 
functions with the spectral function are written by 

Sa{po,p) - r dz P^'-'P}. , (10a) 

Stt 7_oo Pq- z + irj 

Sb{po,p) = r dz , (10b) 

SfiilLOn.p) = r dz^^. (lOc) 

27r y_oo Po- z 

The free spectral function is deriven from the free Dirac equation. Since the Hamiltonian at a zero chemical 
potential is replaced hy H ~ fiN, the free Dirac equation at non-zero chemical potential is written by 

(iTo^o + 7oAi + n • V - m)V-'(.x) = 0. (11) 

Thus, ©, (I5a|) . and (pbl) for the free theory are constructed by V'(a^) satisfying and (ITT]) . S^^^{t,x) obeys 
the equation, 

{ilodo + 7oA* + n • V - m)Sf„^ {t, x) = 0, (12) 
Taking account of ([TTjl . we can obtain the spectral function for the free theory, 

p(°)(p) = 27r((po + m)7o - 7 ■ P + m)e{po + fi)S{{po + fif - |pp - m^). (13) 
where e(po) = 6'(Po) - ^'(-Po)- 



2.2 Formulation (B) 

The definitions of the green functions are the same as ([2]), (l3al) . pbp . and (pc| . However, a fermion field is not 
same. We consider the case that the time evolution is given by 

iP'{t, x) = e'^'i^'it = 0, x)e-'^K (14) 

Then, the KMS condition is given by 

ip'{t,x) = -e^''tlj'{t^il3,x), (15) 

and 

S>'it,x) ^ -e^''S<'it + iP,x). (16) 
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Then, S^'{p) has the relation, 

dt / d?xe'P''S>\t,x) = -e^^P°-''^S<'{p). (17) 

-oo J 

Since the definition for a spectral function is same, S^'{p) and 5<'(p) can be expressed by 

S>'{p) = (1 - npipo - S<'ip) = -nf (po - m)p(p), (18) 

with 

"-f(po - a*) 



An explicit expression for the spectral function is given by [8] 

p'ipo) = 27rZ-i ^ e-«^"-'^^")(l + e^(^"-^'"+^))(m|V^(0)|n)(n|V'(0)|m)<5(po + - i?,n). (19) 
If we perform the Fourier transform, 

POO 

f{po) = / dte^P°'f{t), (20) 

OO 



the retarded green function is 

S'n{p.,p) = r dz P'^'-^Pj. . (21) 
27r J_oo Po- z + iTj 

As can be expected from ^TE\\ and (l?f)) . this does not agree with (jlObp . On the other hand, using the different 
Fourier transform, 



fipo) = / dte^P°*e'^7(0> (22) 

one obtains the result, 

Sr{po,p) - r dz P'^'^P^ . . (23) 

If p{z,p) depends on z + /i, ([221) relates to (I10b|) by the change of variable (see also ([5]) and (fTS)) '). Thus, using 
not ((20)l but (|22p . an equivalent result is derived. This replacement is equal to embedding (|25|) in the Fourier 
transform. 

The free Dirac equation at non-zero chemical potential is 

(i-fodo + ij -V -m)ilj'{x) =0. (24) 

The Dirac equation in the formulation (B) is the same as the zero chemical potential form. From Q, (fTTI) . (|T5ll 
and pi)) . ip{x) and i^'ix) have the relation. 
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Sr'^[{t,x) constructed by ^p'{x) must satisfy the equation, 



(25) 



{ilodo + 7 ■ V - m)S>'^ (t^x)^ 0. (26) 
Taking account of (P^ . we can obtain the spectral function for the free theory, 

p(o)'(p) - 2^(po7o - 7 • P + m)e{po)5ipl n?). (27) 

As pointed out above, this resuh is consistent with ([T5)) . 



3 Cause of the discrepancy 

Obtaining a different resuh is understood from More specificahy, the different results from a plane wave 

solution. Plane wave solutions for (fTT|) and ([Ml) are 

iP{t) ~ e-'^f"-")*, (28) 

and 

^\t) ~ e-'f"*, (29) 

respectively. Thus, there is a difference of ^ for a frequency (energy) between ((TT|) and ([24|). Owing to this, 
after performing the same Fourier transform, the formulation (A) and (B) derive a different result. 
Using (UHl), S>{x) in ^ is written as 



S>{x) = y ^^(z7o(5o - z/i) + Z7 • V + m)e^^*[(l - n;^(£;p - /i))e-^f- - ^^.(ii;^ + ^Ji)e'P% (30) 
where po = Ep = ^Jp^ + m? ■ On the other hand, using Eq. ([2^ . S'>'(x) in ([25]) is written as 

These satisfy each KMS condition. Since there is an extra factor e'^* in (ISO)) , a different result is derived on a 
calculation of a Feynman diagram if one calculates straightforwardly. However, the difference is removed easily 
in the momentum representation using the Fourier transform (1221) . 



4 Example 

4.1 Calculation of the summation by (A) 

We show a procedure of a sum over the Matsubara frequency by using a spectral function[6,. A sum over the 
Matsubara frequency can be done easily by this procedure. We calculate the summation, 
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I = T^uj„iGi3{uj„i,Ei)Gp{un -uj„i,E2), (32) 

m 

where uj„ and LOm are the boson and the fermion Matsubara frequency, respectively. For simpHcity, we use 

Gfj{uJ„i,Ei) ^ — 1 ,„ 2 ' Gp{uJn - UJm,E2) = — ■ 2 ' (^3) 

-(liOm + IJ-r + ^1 -(«W„ - [iLUrn + l^jY + E^ 

(po + ^J■)"fo — ■ p + m in (jl3p is removed. The relation, 

J rdT^^fl^e^"'"-=w„ ^dTG(r,p)e^"'"^ (34) 

can be derived from 



Gp{iuj^,p)= / drG(r,p)e''^"^ 
Jo 

and (HI changed by imaginary time t. Substituting into /, and using ([5]), 



where 



(35) 



M ^ / ^ ; a\'{ko,Ei)a\'{po,E2) 36) 

, 27r 27r za;„ - po - Kq 



(fco, i;i) = 27re(fco + M)<5((fco + m)' " Ef), 
<J^°\po,E2) = 27re(po - m)'5((po - m)' - i^a') 



4.2 Calculation of the summation by (B) 

We consider the same summation ([5^ . However, the same Fourier transform and a spectral function are not 
available. Using the Fourier transform ([2^ . p4p turns into 

dTS{T, i;i)e*("'" = « /'^ ^^ t^g(T,^i) ^,(^^^_,^), ^ ^^^^^^ ^^)g»(-„.-v)r ^ (37) 

Jo Jo Jo 

Substituting this into 
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™ "'0 



, 27r 27r lUn - Po - h 



where 



a(fco,-Bi) = 2^e(fco)(5(fc2 - £;2). (38) 

This resuh agrees with (p6)) . Note that (IT8|) and p7l) are used. Simuhaneous usage of the formulation (A) and 
the modified Fourier transform produces a incorrect resuh. 



4.3 Spectral function in the 1-loop order 

We consider the scalar boson-fermion interaction. We perform the summation without using the Fourier trans- 
form to reduce argument. Ignoring contributions of vertices, the 1-loop fermion self energy in the imaginary 
formalism is given by 

T,{iujn,p) = J j^^Sp{iuji,k)Di3{iuj„ - iuJi,p - k), (39) 

where cij„ and uji are the fermion Matsubara frequency. After replacing 5^ and Dp with (|10c[) . performing the 
sum over the Matsubara frequency by a summation formula without using the Fourier transform, 

f (Pk f°° dzi dz2 1 f , \ , \\ 

J [2t:Y 2tt 2n luj^ - zi ~ Z2 

where nsiz) = l/(e^^ — 1), pp and pB are the fermion and the boson spectral functions, respectively. After 
analytic continuation icOn po + irj to an arbitrary continuous value po, the imaginary part is given by 



1 f d^k f°° 

ImE_R(po,p) = ^ y (27r)3 J dzidz2PF{zi,k)pB{z2,p - k)[nF{zi) + nB{-Z2)]S{po - zi ~ Z2) (41) 

Using free spectral functions (P7)) and Pb\q) — 27re(go)'5(9^ — w^)) 

— TT f d^k 
lmi:n{Po,p) = J (^(^^170 -J-k + m) 

X [(1 - npiEi -p)+ nB{E2))S(po + p - Ei - E2) 
+ {npiEi -p)+ nB{E2))d{po + p - Ei + E2) 

- (1 - npiEi + p) + nB{E2))5{pf) + p + Ei + E2) 

- {npiEi +p)+ nB{E2))S{po + p + Ei - E2)] (42) 
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where El = -^/jfep + m^, E2 = y^\p — fcp + m^. The 1-loop fcrmion spectral function can be obtain by[6][7] 

Pf\p) = -*(S(po + iv) - S(p - ill)) = 2ImSji(p), (43) 

The delta functions in p^p\p) obtained strictly by the imaginary time formalism are different from (50) in Ref 
[5], which used the formulation (B). Thus, the method in Ref [S^ does not correspond to the imaginary time 
formalism. The source of mistake is the Fourier transform (23) in Ref [8^. If one wants to correspond to the 
imaginary time formalism, po must be replaced by pq + /i. However, when one does not consider corresponding 
to the imaginary time formalism, the method in Ref [5] is correct. (As above mentioned, it is nothing more 
than a shift of an energy.) In other words, the difference of /i is absorbed into an external line pq to regard 
^nipojP) as Sfl'(po = Pa + MiP)- In fact, using (0^ . the same retarded green function can be obtained. 



5 Summary 

There are two formulations at non-zero chemical potential. One is the formulation (A) that a Lagrangian 
includes a chemical potential, the other is the formulation (B) that a Lagrangian does not include a chemical 
potential. The former corresponds to the imaginary time formalism more directly. 

It is necessary to note that a calculational procedure in two formulations is different. For example, S^g-^{x) 
in (A) has a factor e*''* as compared with S^^s^{x) in (B). This difference is important when performing a 
summation in imaginary time formalism and calculating a Feynman diagram in the real time formalism, etc. 
A wrong choice makes a mistake. Owing to this, it is necessary to understand the formulation used, combining 
(A) and (B) carelessly is unsafe. In particular, it is important to note that the advanced, the retarded, and the 
thermal green functions expressed by a spectral function are different. 

(B) does not correspond to the imaginary time formalism directly because a shift of an energy exists. This 
shift is removed by modifying the Fourier transform. The modified Fourier transform is equal to embedding 
(1251) in the ordinary Fourier transform. Using the modified Fourier transform, a correspondence between (A) 
and (B) during a calculation becomes easier, (B) agrees with the imaginary time formalism. 

Incidentally, it should be seen that (A) and (B) in the real time functional integral formulation are related 
by the canonical transformJT]. Treating (B) as the canonical transform of (A), a source term has the factor e*^*. 
Thus, the time ordered propagator obtained by a functional derivative in (B) corresponds to the one in (A). 
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